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NONLINEAR AEROELASTIC EQUATIONS OF MOTION OF TWISTED, NONUNIFORM, FLEXIBLE 


HORIZONTAL-AXIS WIND TURBINE BLADES 

by 

Krishna Rao V, Kaza 
The University of Toledo 
Toledo, Ohio 43606 

SUMMARY 

The second-degree nonlinear equations of motion for a flexible, twisted, 
nonuniform, horizontal-axis wind turbine blade are developed using Hamilton^ 
principle* The derivation of the equations has its basis in the geometric 
nonlinear theory of elasticity, and the final equations are consistent with 
the small deformation approximation in which the elongations and shears are 
negligible compared to unity and the square of the derivative of the 
extensional deformation of the elastic axis is negligible compared to the 
squares of the bending slopes* A mathematical ordering scheme which is 
consistent with the assumption of a slender beam is used to discard some 
higher-order elastic and inertial terms in the second-degree nonlinear 
equations* The blade aerodynamic loading which is en$>loyed accounts for both 
wind shear and tower shadow and is obtained from strip theory based on a 
quasi-steacty approximation of two-dimensional, incompressible, unsteady, 
airfoil theory* The resulting equations have periodic coefficients and are 
suitable for determining the aeroelastic stability and response of large 
horizontal-axis wind turbine blades. 
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INTRODUCTION 


The recently renewed efforts in wind power are due to its prospective 
uses as an alternative energy source. As a result of these efforts, several 
wind turbine projects have been initiated by NASA Lewis Research Center as a 
part of the Department of Energy's (DOE) overall wind energy program. To make 
wind energy cost effective, wind turbines substantially larger than the 
existing lOOkW Mod-0 which has a rotor diameter of 38 meters are being 
studied. However, as the rotor diameter increases, blade flexibility and 
hence susceptibility to aeroelastic instability also increase. Furthermore, 
efficient construction and operation of wind turbines require that the 
vibratory loads and stresses on the rotor as well as on the combined 
rotor-tower system be reduced to the lowest possible levels. Thus, 
aeroelastic and structural dynamic considerations have a direct bearing on the 
manufacture, life, and operation of these large wind turbine systems. 

Although the structural ctynamic and aeroelastic technology used to develop 
helicopter rotors appear to be adequate for the development of wind turbine 
machines, this technology has to be transformed from helicopter rotor 
applications to wind power applications, and additional studies have to be 
conducted to determine the effects of the parameters peculiar to wind power 
machines on the aeroelastic and structural dynamic behavior. 

Several aeroelastic considerations are common to both the wind turbine 
and helicopter blades. These include flap-lag-torsion, flap-torsion and 
flap-lag instabilities, stall flutter, and torsional divergence. The wind 
velocity gradient due to the Earth's boundary layer and gravity loads in the 
case of a wind turbine rotor, and forward velocity in the case of a helicopter 
rotor lead to timewise periodic coefficients in the equations of motion. 
Several previous studies have considered the helicopter blade and developed 
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the nonlinear aeroelastic equations of motion. Recently, reference 1, though 
primarily a study of the flap-lag <fynamics of rigid articulated helicopter 
rotor blades, contains a cursory examination of the elastic blade. References 
2 and 3 more completely examined the basis of the nonlinear aeroelastic 
equations. In particular, attention was directed at establishing the 
expressions for the nonlinear curvatures and the nonlinear transformation 
matrix between the undeformed and deformed blade coordinates. The resulting 
equations were compared with some of those in the literature. These 
comparisons indicated several descrepancies with the results of reference 3, 
particularly in the nonlinear terms. The reasons for these descrepancies were 
explained in reference 3. 

For wind turbine blades, reference 4 presented a set of nonlinear 
equations of motion. An examination of these equations reveals that reference 
4 fails to obtain several nonlinear elastic and aerodynamic terms which are of 
the same order as those retained. It appears that these terms were not 
obtained in reference 4 for two reasons: (1) an incorrect torsional curvature 

egression and (2) linearizing the resultant transformation matrix between the 
undeformed and deformed blade coordinates while developing nonlinear equations 
of motion. 

Based on the considerations and discussions in references 1 and 2, 
reference 5 apparently sought to redevelop the nonlinear equations of 
equilibrium for rotor blades. However, the resulting equations are missing 
several elastic terms which were well established in the literature and some 
nonlinear terms as well. It is interesting that several third-degree 
nonlinear terms were retained. These general equilibrium equations were used 
to study aeroelastic stability of a single wind turbine blade in reference 6 
and a coupled rotor/support system in reference 7. Nonlinear aeroelastic 
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equations of motion of a single wind turbine blade were also developed in 
reference 8 using the Newtonian method. Several nonlinear elastic and 
aerodynamic terms are missing in this reference because of the use of an 
incorrect expression for the torsional curvature and a partial linearization 
of the resultant rotational transformation between the coordinates of the 
undeformed and deformed blade. In view of the differences in the equations 
existing in the literature, it is felt that a comprehensive development of the 
nonlinear aeroelastic equations of motion of an elastic horizontal-axis wind 
turbine blade is required. The basic ingredients of such a development were 
presented in reference 9. This report documents the details of the 
development. 

The derivation of the nonlinear equations of motion herein follows 
essentially along the lines of reference 3. In this reference, the pretwist 
together with control inputs of the blades were combined with the elastic 
twist for simplicity, following a common practice in the helicopter blade 
literature. Physically, the pretwist is present in the blade even before the 
deformations. Thus, this report will include the rotation due to pretwist 
with the control inputs and inqpose this rotation first while inposing the 
rotational transformations. 
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SYMBOLS 


a airfoil lift-curve-slope 

A cross-sectional area 

A U ,A V ,A W generalized aerodynamic forces per unit length in X, % Y, Z 

directions , respective ly 

A^ generalized aerodynamic moment per unit length about 

elastic axis 

b number of blades; also blockage factor to account for 

tower shadow 

B V/ B T# B a ,B5 W ,B a boundary terms arising from strain energy, kinetic energy, 
y D 

gravity, material damping, and aerodynamic forces, 
respectively 
section constants 
blade chord 

airfoil profile drag coefficient 


B 1 9^2 

c 

c d 

o 


C(k) 


c 1» c 2 


d^ (i=1 , 2 , • • • 5) 


D 

e 


e A 


E 

E* 


Theodorsen's circulation function 
aerodynamic constants 

notation used in writing the virtual work associated with 
material damping in concise form 
airfoil profile drag per unit length 
chordwise offset of mass centroid from elastic axis 
(positive when in front of elastic axis) 
chordwise distance of area centroid of cross section from 
elastic axis (positive when in front of elastic axis) 
Young's modulus 

coefficient of internal friction in tension 
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e v e y 3 /ez 3 


e X' e Y' e Z 


unit vectors along x 3 , y 3, Z3 axes, respectively 
unit vectors along X, Y, Z axes, respectively 


eX l' eY l' eZ I 


unit vectors along X-j, Y-), Z-| axes, respectively 


e x» e n' e £ 

F x 3 ,p y 3 /Fz 3 


F x' F n' F c 


F gu /F gv ,F gw 


F g<t> 

g 

^1 • ^2 • * • • ^6 
G 
G* 
h 


H 




j 


unit vectors along x, r\, C axes, respectively 
components of aerodynamic forces per unit length in X3, 
y3, Z3 directions, respectively 
components of aerodynamic forces per unit length in x, 
h, C directions, respectively 
corqponents of gravitational forces per unit length in X, 
Y, Z directions, respectively 
gravitational moment per unit length about elastic axis 
gravitational constant 

notation used in writing the gravitational forces 
shear modulus 

coefficient of internal friction in shear 
vertical velocity of two-dimensional section normal to 
f ree-stream 

height of rotor hub center above ground 
generalized inertia forces per unit length in X, Y, Z 
directions , respective ly 

generalized inertia moment per unit length about elastic 
axis 

area moments of inertia about r\ and ? axes, respectively 
torsional section constant 


k 


reduced frequency 



k A 

^( 1 = 1 , 2 , ... 6 ) 


l^m^n* 

L 




m 

<J> 

R 

*1 

*o 

r i 

r 0 

s 

S 1 

Si (i=1,2,,..10) 


T 

T A 

T C 


polar radius of gyration of cross-sectional area about 
elastic axis 

notation used in writing the variation of kinetic energy 

polar radius of gyration of cross-sectional mass about 
2 2 2 

elastic axis (1^ = k^ + k m ^) 

direction cosines (i - 1, 2 , 3) 
aerocfynamic lift per unit length 

aerodynamic pitching moment per unit length about deformed 
elastic axis 

mass of blade per unit length 
generalized dartping forces 
length of the blade 

position vector of a point on the elastic axis after 
deformation 

position vector of a point on the elastic axis before 
deformation 

position vector of a point after deformation 
position vector of a point before deformation 
coordinate along the undeformed elastic axis 
coordinate along the deformed elastic axis 
notation used in writing the variation of the strain 
energy 

kinetic energy 

blade tension due to aero<fynamic forces 
blade tension due to centrifugal forces 
blade tension due to gravity 
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U R ,U T ,Up 


u 

Up 


u,v,w 


Ug/ Vqi w g 

V 

V H 

V m Ul> 

Vl Z I 

Vx 3 y 3 Z 3 

^xn c 


V i 

w 

W A 

W D 

w g 

XYZ 


radial, tangential, and perpendicular components of 
velocity for blade airfoil section 
resultant of U T and Up 
radial foreshortening of elastic axis 
deformation of elastic axis in X, Y, and Z directions, 
respectively 

components of gust velocity in Xj, Yj, Zj directions 
respectively 
strain energy 

mean wind velocity at height H 
mean wind velocity at xj 

aerodynamic velocity vector e^qpressed in XfYjZi coordinate 
axis system 

relative velocity of point on elastic axis expressed in 
X 3Y3 Z 3 coordinate system 

relative velocity of point on elastic axis expressed in 
xtiC coordinate system 
induced velocity 

work done by aerodynamic, structural damping, and 
gravitational forces 
work done by aerocfynamic loading 
work done by structural damping 
work done by gravitational forces 

coordinate system with the origin at the hub center 
which rotates with blade such that X-axis lies along 
initial undeforraed position of the blade elastic 
axis 
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Xl*lZi 

XftYazn 

xn? 

x oyo z o 

X 1/Y1/ Z 1 

X 3Y3 Z 3 

[T] 

[Tr] 

C^ij] 

a 

&pc 

3 i, Cl# e 1 


^xx' ^xn' 


inertial axis system with the origin at hub centerline and 
Zj-axis normal to the hub plane 
hub-fixed axis system rotating about Zj-axis with 
angular velocity ft 

principal axis system obtained by rotating about xg-axis 
with an angle 0 

blade-fixed axis system at arbitrary point on elastic axis 
before deformation 

coordinates of a point (which was at xnC in the un deformed 
blade) in the deformed blade 

blade-fixed orthogonal axis system in deformed configura- 
tion obtained by rotating xnC; X3-axis is tangent to the 
deformed elastic axis 

transformation matrix between xnC and X3y3Z3 coordinate 
systems 

transformation matrix between xn? and XjYjZj coordinate 
systems 

Green 1 s strain tensor 

airfoil section angle of attack, a = tan"^Up/U»p 
notation used in writing the derivative of the 
displacement of a point 
angle of built-in coning (precone angle) 

Eulerian-type rotation angles between xti? and 
X 3Y3Z3 

engineering strain conponents 


6 ( ) 


variation of ( ) 
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e xx' e xn* e x? 
n 


e 

x 
e 

0 c 

0 pt 

x i 

V 

X H 

\iq* ^vg' \rq 
P 

°xx' a xn' a x£ 
♦ 

+ 

♦o 

0) X ,0>y ,0)2 

3 y 3 3 


small parameter of the order of the bending slopes; the 
airfoil section pitch angle with respect to free-stream 
velocity; extensional cotqponent of Green's strain tensor 
on the elastic axis 
tensor strain corrponents 

sectional coordinate along major principal axis; also 
constant for mean velocity calculations 
sectional coordinate normal to n-axis 
non dimensional coordinate along blade axis, x = x/R 
total geometric pitch angle, 0 = 0 pt + 0 C 
collective pitch 

built-in twist (pretwist), positive when leading edge is 
upward 

induced flow ratio, X^ftR = v^ 
inflow ratio, V* R = V n (x!) 
inflow ratio, XflftR = V H 

inflow ratios, X ug ftR = u g sin u> g t; X vg ftR = v g sin u) g t; 

X wg ftR = w g sin o) g t 

mass density of the blade; also mass density of air 
engineering stresses 

angle of twisting deformation about elastic axis, 
positive when leading edge is upward 
blade azimuth angle measured from downward position in the 
direction of rotation 

angle between mean wind velocity V m and rotor axis Zj 
torsional curvature (total rotation rate about X 3 -axis) 
and bending curvatures, respectively 
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curvature vector before deformation 


“xnc 


U). 


* 3 * 3*3 


Q 

( >c 

( >NC 
[ ] T 
( ) 

( )' 


curvature vector after deformation 

angular velocity of xtiC coordinate system 

rotational speed of the wind turbine rotor 

circulatory aerodynamic term 

noncirculatory aerodynamic term 

denotes transpose of matrix 
3 

time derivative, -- ( ) 

3t 

3 3 

space derivative, — ( ) = — ( ) 

3s 3x 
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MATHEMATICAL MODEL AND COORDINATE SYSTEMS 


The mathematical model chosen to represent the wind turbine blade 
consists of a straight, slender, variably twisted, nonuniform elastic blade. 
The elastic axis, the mass axis, and the tension axis are taken to be 
noncoincident; the elastic axis and the feathering axis are assumed coincident 
with the quarter chord of the blade. The effect of warping is assumed to be 
small and is neglected. The generalized aerodynamic forces account for tower 
shadow, wind shear, and gusts, and are calculated from strip theory based on a 
quasi-stea<fy approximation of two-dimensional, inconpressible, unsteady, 
airfoil theory. 

Several orthogonal coordinate systems will be enployed in the derivation 
of the equations of motion; those which are common to both the dynamic and 
aero<fynamic aspects of the derivations are shown in figures 1 to 3. The axis 
system XjYjZj shown in figure 1(a), is fixed in an inertial frame with the 
origin at the center of the hub. The axis system XYZ, shown in figure 1(b), 
is obtained by rotating XjYjZj system about Zj— axis by an angle i|> = ftt and 
then XftYftZjj system about the negative Yjj-axis by an angle 0 pc , the angle of 
built-in coning angle. The point on the cross section through which elastic 
axis passes is given by the intersection of the Y and Z axes. Let Xoy 0 z 0 be 
axes fixed to the blade at an arbitrary point on the elastic axis of the blade 
so that before deformation XqYqZq are parallel to XYZ respectively. The ti and 
£ axes (fig. 2) with the origin at the elastic axis of the cross section are 
principal axes and are inclined to the Y and Z axes by an amount 0. The 
geometric pitch angle is given by 6 = 0 pt + '.0 C where 0 pt is the built-in twist 
angle (pretwist) and 0 C is the collective pitch angle. 

The variables defining the configuration of the deformed blade are shown 
in figure 3. When the blade deforms, the elastic center of an arbitrary 
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section deforms an amount u in the X direction, v in the Y direction, and w in 
the Z direction. The section rotates about the principal axes due to bending 
in addition to twisting an amount <|> about the elastic axis. The final 
position of the triad xn5 after deformation is denoted by X 3 y 3 Z 3 » 

HAMILTON'S PRINCIPLE 


The equations of motion will be derived using the extended Hamilton's 
principle as in reference 3 in the form 


r t. 

} 1 (6T - 6V + SW) dt = 0 

to 

where 


(1) 


<5W = 6W D + 6w g + 6 W a (2) 

In equation (1), T is the kinetic energy, V is the strain energy, and 6W is 
the virtual work done by the damping, gravitational, and aerodynamic forces. 
In the following sections ejqplicit expressions for 6T, 6V, and 6W in terms of 
the variables u, v, w, and <J> and the blade sectional properties will be 
developed. 


STRAIN ENERGY 


The expression for the strain energy of the blade in terms of stresses 
and engineering strains is 

1 ,R , 

v = “ Jq J a ( a xx^xx + a xrpfxn + °x?Yx 5) dx (3) 

where, using Hooke's law. 


*The coordinates s and x are used interchangeably. 
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°XX ™ E ^xx 

°xn = G Yxn (4) 

a xc = G YxC 

Assuming small strains, the engineering strains are related to the conponents 

of the strain tensor according to 


Yxx - 

e xx 


^xn = 

2e xn 

(5) 

Yx? = 

2 e x? 



To develop the explicit expressions for strains, the expressions for the 

curvatures u>„ , u> v , u> z and the transformation matrix [T] between the 
^3 y 3 3 

coordinate systems xn5 and X 3 y 3 Z 3 are required in terms of the variables u, 
v, w, <j> and the blade sectional properties. These expressions are developed 
in Appendix A following the procedure described in reference 2 and are given 
by 

,2 ,2 

v' w' 

w x ™ <J>' + 9' (1 - ) - (v' cos 9 + w' sin 0) 

3 pt 2 2 

* (-v' sin 0 + w' cos 0)' 


= -w' ' (cos 0 - (j> sin 0) + v' '(sin 0 + <f> cos 0) (6) 


w z = v' ' (cos 0 - <j> sin 0) +w ,, (sin 0 + <j> cos 0) 
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[T] = 


'2 '2 

V * W * 

2 2 


v'cos 0 + w' sin 9 | -v' sin 9 + w' cos 9 


I 


-v' (cos 9 - <(> sin 0) 


-w' (sin 0 + <|> cos 0) 


v' (sin 0 + <J> cos 0) 


-w' (cos 0 - <j> sin 0) 


1 - (v'cos0+w'sin0) 2 
2 

4> 2 


-♦ 


$ - (v' cos0+w'sin9 ) 


• (-v' sin0+w'cos0 ) 


1 - 


r v 


(-v'sin0+w'cos0 ) 2 


(7) 


Using equations (6) and (7), the second-degree nonlinear expressions for 
the strains are developed in Appendix B, and are 

<J > 2 

= e vv = u' + (n 2 + 5 2 ) ( — + 9' $' ) - (v' ' + w' '<|))(bcos 0-5 sin 0) 


'XX XX 


2 pt 

- (w' ' - v* ' <j>) (n sin 0 + 5 cos 0 ) 


(8a) 


y xr , = 2e xn = -5 [<f>' - (v' cos 0 + w' sin 0) (-v' sin 0 + w' cos 0)' 


Yx 5 = 2e x5 


n ir 


v'2 w'2 

-0* ( + )] (8b) 

pt 2 2 


- (v' cos 0 + w' sin 0) (-v' sin 0 + w' cos 0)' 


v' 2 w' 2 

-0' ( + )] (8c) 

pt 2 2 

In the above expressions/ several higher-order terms have been discarded based 
either on considerations related to the small deformations Level I 
approximation/ as discussed in reference 2, or on considerations related to 
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the approximations which can be made because of the assumed slenderness of the 
blade, as discussed in Appendix C. Formal retention of higher-order terms in 
the expression for strain components is not a problem. However, these 
higher-order terms will lead to higher-order elastic terms in the final 
equations of motion. Thus, discarding these higher-order terms in the 
expressions using the considerations of Appendix C simplifies the subsequent 
algebraic manipulations. 

Taking the first variation of V as given in equation (3) and using 
equation (4), yields 


6V = J E // YxxSYxxdn dQ dx 
0 A 

+ /*<;// ( Yxn fi Yxn + Yx^Yx*;) dT1 dx (9) 

0 A 

Substituting equation (8) into equation (9), taking the variations, and 
integrating over the cross section leads to 

R 

6V - / (s 1 6u * + s 2 6$ + s^ 1 + s 4 6v' ' + s_6w' ' + s 6<j>» + s 6v* 

0 7 

+ Sg6v' ' + s 9 6w' + s 1() 6w") dx (1( 

where 


= EA [u* + k 2 ( d>» 0 • + - <b * 2 ) 

A pt 2 

- e A (v’ ' + ij)w" ) cos 0 + e A (<j>v' • - w" ) sin 6] 

= EAe^u' (v* 1 sin 0 - w 1 1 cos 0) 

+ EB 2 ( v 1 1 0 1 4> f sin 0 - w' , 0* <j)' cos 0) + v' 'w 1 ' (EI CC - EI nn )cos 20 

pt pt ,|M 

+ w" 2 (EIjr 5 - EI nn )sin 6 cos 0 + v’ ' 2 (EI nn - EI ?? ) sin 0 cos 0 
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s 3 = EAk 2 u’(0' + <J>') + EBt [6’ <J>' 2 + 9’ (<J>'0 ' + - <b’ 2 )] 

A pt pt pt pt 2 

+ EBo [0 1 ((bv' ' - w' ' )sin 0 - (j>'v' ' cos 0 - (b’w' ' sin 8 
pt 

-0' (v' ' + <j>w' ' ) cos 0] 

pt; 

8 4 = EAeRu'(<J> sin 0 - cos 0) - EB2<f>'0' cos 0 

pt 

+ v* ' [El-pj^tsin 2 © + <J) sin 20) + EI^ (cos 2 0 ~ <J> sin 20)] 

+ w’’[(El££ -EI^rjMsin 0 cos 0 + <J> cos 20)] 

s 5 = -EAeRu'(<(> cos 0 + sin 0) - EB2<b , 0' sin 0 

pt 

+ v''[(El££ - EI nT1 )sin 0 cos 0 + (j>(EI^^ - EI r ^ r( )cos 20] 

+ w" [EI nn cos 2 8 + EI 55 sin 2 6 + - Eltsin 20] 

s 5 - GJ[<J>* + (v'v* 1 - w’w* 1 ) cos 0 sin 0 - vV 1 cos 2 0 

v' 2 w' 2 

+ w'v* 1 sin 2 0 + 0* (- — cos 20 cos 20 + v’w’ sin 20)] 

pt 2 2 

s 7 = GJd) 1 (v' 1 cos 0 sin 0 - w' ' cos 2 0 + v’0 1 cos 20 + w 1 © 1 sin 20) 

pt pt 

8 8 = GJ^'Cv* cos 0 sin 0 + w* sin 2 0) 

s q = GJ6' (-w 1 1 sin 0 cos 0 + v 1 1 sin 2 0 + v* 0 1 sin 20 - 0 1 w’ cos 20 ) 

pt pt 

s 10 « -GJ^'tw' cos 0 sin 0 + v* cos 2 0) 

The sectional properties appearing in equation (11) are defined as follows: 

A ® // dn dC Ae = / / ri dri dc 

A 

inn = // s 2 dn Hi = // n 2 dn dc 

Ak 2 = // (n 2 + ? 2 ) dn dc j = // (n 2 + i 2 ) dn <15 (12) 

A 

2 

B i “ // (n 2 + l 2 ) dn d^ b 2 « // n(n 2 + l 2 ) dn d^ 
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Since the cross section is assumed symmetrical about the n axis, the 
following integrals are zero: 

// 5 dn d? = 0 // nc dn dc = o 

//■ c (n 2 + C 2 ) dn dc = o (13) 

Integrating equation (10) by parts, the resulting expressions can be put 
in the form 


6V = / (S u Su + S v 6v + S w 6w + S* 6<|>) dx + By 
0 


(14) 


where the generalized elastic forces S^, S v , S,,, and S<j,, to second-degree, are 
given by 

S,, = -s’ 


s' ' - s 


' + s' ' 
7 8 


(15) 


J w 


,i i - 


s' + s ' ' 
9 10 


Srf, = s - s' - s' 
9 2 3 6 


and the boundary term B v is given by 

B „ = rs 6u + (s + s )6v' + (- s * + s _ s' )6v 
V 1 4 8 4 7 8 

+ (s + s )6w' + (-s' + s - s' )6w + (s + s > 6 4>1 J (16) 
5 10 5 9 10 3 6 0 


KINETIC ENERGY 


The expression for kinetic energy T is given by 

1 R . . dri dr 1 

T = - III p . dr) d£ dx 

2 o A dt dt 


( 17 ) 
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and its variation, integrated between t 0 and t, , is given by 


( 18 ) 


tl t, » dr, dr, 

f St = / f II P • 5 dn d£ dx dt 

to to o A dt dt 

In the above equation, the absolute velocity of the mass point is dr,/dt and 

is defined by 


dr, - - - 

= r, + to x r, (19) 

dt 


where to is the angular velocity of the xn5 coordinate system and r, is the 
position vector of the mass point expressed in terms of the unit vectors of 
the xnC system. The angular velocity u) is obtained by projecting £) along the 
x, tl/ and 5 directions and is given by 


to = npp C e x + & sin0 e,, + ft cos0 e^ 


( 20 ) 


In the above equation, the precone angle 0p C is assumed to be small. The 
position vector r, is given by 


+ 


= xe x + (u - Up)e x 

(v cos 0 + w sin 0 )e^| + (-v sin 0 + w cos 0)e^ + [T] 


n >(21) 


Substituting equation (A38) into equation (21), the conponents of the position 
vector to second-degree are given by 


c , = x + u - Up 

- (v’ + w'ij))(n cos 0 - ? sin 0) - (w' - v'<J>)(h sin 9 + C cos 0) 

(v' cos 0 + w' sin 0)2 <|>2 


y, = v cos 0 + w sin 0 + [1 - 


■] Tl - C<j> 
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— v sin 0 + w cos 6 


(22) 


'1 = 


+ [<t> 

+ [1 


(v' cos 0 + w' sin 0)(-v' sin 0 + w' cos 0)] n 
(-v' sin 0 + w' cos 0)2 <J> 2 


Differentiating r -j with respect to time according to equation (19), the 
absolute velocity of the mass point is 
dr-i 

_ (x* + zjl sin 0 - y^ cos0 fi)e x 

dt X 

+ (Y., " z ! n 3 pc + x ! cos6 fi) ® n + (S i + “ x 1 ft sin0)e C (23) 


Substituting equations (22) into (23) and the result into equation (18), 
integrating by parts over time where necessary, and then integrating over the 
cross-section, the variation of T can be put in the form 

5T = /* <V U " k 1 6U F 
0 

+ k 2 <Sv' + k.j6v ” k^Sw' - k,.Sw + k 6 6<|>) dx (24) 

where, consistent with the ordering scheme given in Appendix C, are 
#* 

k 1 = - mu + 2mft(v - e$ sin 0) 

+ mQ 2 (x + u - U F - ev' cos 0 - ew' sin 0) 

• • •• 

- mft 2 3p C (w + e sin 0 + e<j> cos 0) + me (v' cos 0 + w' sin 0) 

k 2 = mft 2 e<|>x sin 0 - 2mefiv cos 0 - meft 2 x cos 0 

• • •• 

k 3 = m« 2 (v + e cos 0 - e<J> sin 0) - mv + me <)> sin 0 + 2mft3p C w 

- 2mft(u - Up - ev* cos 0 - ew 1 sin 0) 

« mft^e<f>x cos 0 + 2mftev sin 0 + ex sin 0 
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= mw + me <f> cos 0 + 2mftgp C v + m^ 2 3p C x 

^6 ~ ‘ -w ' (2mftev cos 6 + mfi 2 ex cos 0) - mQ 2 e$v cos 0 

- mft 2 <|>(k m - k 2 ) cos 2 0 - mft 2 (k 2 - k 2 ) sin 0 cos 0 
2 m-| m 2 m^ 

+ me<f>v cos 0 + 2mftev , v sin 0 + mft 2 exv' sin 0 
-mfl 2 ev sin 0 + me<J>w sin 0 + raft 2 3p C e<J>x sin 0 

+ 2m ft [e sin 0 u - U F ) - (k 2 - k 2 )v' sin 0 cos 0 

m 2 


•• 

- w’ (k 2 sin 2 0 + k 2 cos 2 0)] + mev sin 0 
m 2 m 1 



ex cos 0 


- mew cos 0 - 2mft3p C ev cos 0 


The sectional properties appearing in equation (25) are defined as follows: 


m 



P dn d£ 


me 



pTl dn dc 


mk 2 = // PC 2 dn d £ mk 2 = // pn 2 dn d£ 

m-j A m 2 A 


(26) 


k 2 = k 2 + k 2 
m m 2 


From symmetry of the cross-section about the n axis, the following integrals 
have been set to zero: 


// PC dn dc - 0 


// pnc dn dc = o 


(27) 
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Since Up is a function of v' and w', the term involving SUp in equation (24) 
requires separate treatment. Using equation (A16), the second term in 
equation (24) can be written in the ejqpanded form 


= / R ^1 ^ w '^ w ' + v '^ v ' ) dx] dx 

0 0 0 

which can be further rewritten as 

/ R k-|^ U p dx = J R ( J R dx) (w'fiw' + v'6v' ) dx 
0 Ox 

Defining the tension T c as 



x 


equation (29) can be written as 

/ X dx = / R T c (w'6w' + v'6v' ) dx 

0 0 


(28) 


(29) 


(30) 


(31) 


Integrating equation (24) by parts, the resulting expression can be put in the 
form 


6T = J* (I u 6u + I v 6v + 
0 

where the generalized inertia 


ft 

H 

K 1 = -T ' 
C 

T a 

A V 

-k' + k + (T c v' ) 
2 3 

II 

* 

H 

k' “ k + (T c w' ) ' 
4 5 

i<j> = 

k 6 


and the boundary term B x by 


I w 6w + 1^6$) dx + B x (32) 

forces I u , I v , I w , and are given by 

1 (33) 
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R 


Bip 


(k 2 - T c v' )6v - (k 4 + T c w' )Sw 

* 0 


(34) 


VIRTUAL WORK DUE TO MATERIAL DAMPING 


The virtual work due to the dissipative forces associated with structural 
(material) damping can be ejqpressed in the form 

6W D = ^ Q Dk 6q k 05) 

where Q^is generalized dancing force associated with the k^* 1 dependent 

variable and is the variation of the k^h dependent variable* In the 
present development the generalized dairying forces accounting for the 
dissipation of energy due to material dancing will be taken to be those 
consistent with the assumption of a material which exhibits a linear 
viscoelastic behavior* This theory assumes that the stresses are linear 
functions of the strains and strain rates. Such a behavior is analogous to a 
spring and a dashpot in parallel, and a model which exhibits such a behavior 
is often termed a Kelvin-Voigt solid in the literature. This model was used 
in reference 3* For the stresses and strains of interest herein, these 
constitutive relations have the form 


T 

XX 

- E ?xx 

+ 

E *Yxx 


T xn 

■ G Yxn 

+ 

G *Y xn 

(36) 

T x? 

= G Yx? 

+ 

G *Y X £ 



where E and G are Young's modulus and the shear modulus, respectively, and E* 
and G* are coefficients which take into account internal damping of the 
material in tension and shear, respectively. The first term on the right hand 
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side of each of equations (36) contributes to the usual elastic strain energy 
and have alreacty been treated in an earlier section. Considering only the 
dissipative terms in equations (36), the virtual work of the structural 
dissipative forces can be written as 


,R , , » 

fiW D " ~ J E * I J Yxx^xx dT1 631 

0 A 

Re • 

- J «* ff (y xn^xj) + WW d * d c < 37 > 

0 A 

The result given in equation (37) is general. However, because of the lack of 
knowledge as to the distribution of dancing, only the direct danping terms are 
generally retained in practice. Thus, off-diagonal terms accounting for 
danping coupling between the dependent variables which arise from equation 
(37) are taken to be zero and only the direct danping terms associated with 
the dependent variables are retained. In addition to adopting this expedient 
in the present development, it will also be assumed that a first approximation 
to the direct damping terms can be obtained by retaining only the linear 
damping terms in the final equations of motion. Thus, it is sufficient to 
retain terms up to only the first-degree in the expressions for the strains. 

To first degree the resulting strain expressions are 

Y xx 3 u ' + + 5^) 9 ' - v* ' (H cos 6 - 5 sin 6 ) - w' ' (h sin 0 + £ cos 0 ) 

pt 

Y xn = ?♦' < 38 > 

YxC 3 

Substituting equations (38) into equation ,(37) , integrating over the cross 
section, and retaining only the linear direct damping terms leads to 


SW D 



0 


(d-| 6 u' + d 2 < $ 4 ), + d 3 <Sw n + d 56 <j>')dx 


(39) 
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where 


d -j = E*Au' 

d 2 = E*Bi 6 1 2 
pt 

d 3 = E*(I^^ cos 2 0 + 1^^ sin 2 0)v' ' (40) 

d 4 *» E*(I^^ sin 2 0 + 1^ cos 2 0)w' ' 

d 5 = G*J4>' 

Integrating equation (39) by parts, the generalized damping forces Qp , Qd • 

u v 

Q d , and Qq become 
w <f> 


Q d = d' = (E*Au ' ) ' 
u 1 

Q d » -d^' = - [EMI 55 cos 2 0 + I nT1 sin 2 0 )v' ' ]" 

Q d = -d" = - [E*(Irr sin 2 0 + I nr , cos 2 0)w ,, ] ,, (41) 

w 4 

Q d » d' + d' = (E^B-j 0 * 2 <(> * ) ' + (G*j5>')' 

4> 2 5 pt 

and the boundary term B$ w becomes 

D 


b 6 w d - - d i 6u l +< v 6v 


R R |R 

- d 3 «Sv' + d 4 '6w 

0 0 '0 


) R i R 

- (d 2 + d^ ) 6 <J) 

0 0 


VIRTUAL WORK DUE TO GRAVITY 


(42) 


The virtual work due to the gravitational forces of the blade can be 
expressed in the form 


<SW g = / R // A pg . 6r 1 dn d ? dx 


(43) 
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where g is the gravitational acceleration vector. It is obtained by 
projecting the gravity, g, which acts along the Xj axis as shewn in figure 
1(b), along the x, n, and C, and is given by 

g = g [cos tex - (cos 0 sin i|i + 0p C sin 0 cos lJOe^ 

+ (sin 0 sin <|i - 0p C cos 0 cos 'JOe^] (44) 

In the above equation the precone angle 0p C is assumed small and hence cos 0p C 
is replaced by one and sin Ppc b y Bpc* 

The position vector, r-j of a point on the blade after deformation is 
given by equation (21). Taking the virtual variation of equation (21), 
substituting the result into equation (43), and integrating the result by 
parts, one obtains 

, R 

= / (Fg^u + FgvSv + FgvSw + F g< j,6<J>) dx + Bg (45) 

0 

= g 

i 

■ V <T s v ' 1 ■ < 

= g + (TqW* ) - g' 

3 6 

= g 

4 

" f" 9 , 4 * 

X 1 

= mg cos 
- -mg sin ij> 

= -mg3 pc cos 


5w„ 


where 


gu 


gv 


F, 


gw 


’g<t> 
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^4 “ mge cos*(v* sin 0 - w' cos 0) + mge<J>(cos 0 sin 

+ (5p C sin 0 cos if) + mge(sin 0 sin iJj - f5p C cos 0 cos Tp ) 
g 5 «* mge cos \|>(-cos 0 + $ sin 0) + mge(cos 0 sin iji 

+0p C sin 0 sin \|» ) . (v' cos 2 0 + w' sin 0 cos 0) 

9 q = mge cos ij)(— 4* cos ® + sin + nige(cos 0 sin ip 

+ 3p C sin 0 cos ip) . (w' sin 2 0 + v* sin 0 cos 0) (46) 

The boundary terms B g are given by 


= ( g 5 ** TgV' )6vl 


+ ( g 6 - TgW' )6w 


(47) 


It should be pointed out that all the terms of order 0(e 3 ) and higher are 
neglected in the expressions for gravitational forces. 


GENERALIZED AERODYNAMIC FORCES 

The aerodynamic forces will be generated from two-dimensional, 
incompressible, quasi-steacfy, strip theory in which only the velocity 
components perpendicular to the spanwise axis of the deformed blade (the X3- 
axis) are assumed to influence the aerodynamic loading. In calculating the 
velocity components, the effects of wind shear and tower shadow as well as 
gusts are included. Account is taken of the pulsating free-stream velocity 
V(t) associated with a rotating blade enploying Greenberg's extension of 
Theodorsen's unsteady theory (ref. 10) for determining the aerodynamic lift 
and pitching moment acting on the blade. The resulting expressions are 
specialized to the case of quasi-steacfy flow by setting Theodorsen's 
circulation function to unity. The classical blade element momentum theory 
can be used to calculate the steady flow induced by the rotor. 

In the present application of Greenberg's theory, the airfoil is taken to 
be pivoted in pitch about the ^ergdynamiq center at the quarter chord and to 
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be executing harmonic motions in pitch (e(t)) and plunge (h(t)) while Immersed 
in a pulsating airstream V(t), as shown in figure 4. The lift and moment 
acting on an elemental section of the blade may be expressed in terms of the 
circulatory and noncirculatory components as 

L = + LflQ (48) 

M = Me + %c 


Assuming that the blade elastic axis is coincident with the aerodynamic center 
at the quarter chord, the individual conponents of equation (48) follow from 
reference 10 and can be written as 


0 •• •• 

1 c^ • c 

Ljjc = - pa — (h + ve + Ve + - e) 
'24 4 


L C 


1 . c • 

= - pacV(h + Ve + - e) 

2 2 


(49a) 


(49b) 


m NC 

MC 


1 c 2 . 3c ~ 

- pac (-) (Ve + h + — e) 

2 4 8 


1 c 2 
pac (-) 2Ve 

2 4 


(49c) 

(49d) 


In the course of arriving at the circulatory terms in equation (49), the 
quasi-steady approximation has been introduced by setting the reduced 
frequency k to zero, in consequence of which Theodorsen's circulation function 
C(k) assumes the value of unity. The noncirculatory lift and moment are 
associated with apparent mass forces and are oftentimes discarded in rotor 
blade applications* Note that Greenberg's modification (i.e., a pulsating 
stream in which V * 0) appears only in the noncirculatory expressions for the 
lift moment* Hence/ if one assumes, a priori, that apparent mass forces will 
be neglected, there is no Greenberg 1 2 s modification* 
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The lifts and moments given in equation (49) must now be expressed in 
terms of U R , U T , and Up, the radial, tangential, and perpendicular velocity 
conponents relative to a point on the elastic axis of the airfoil (fig. 5). 
Now the expression in parentheses of equation (49a) for is the downward 
acceleration of the mid-chord point of the airfoil, and the expression in the 
parentheses of equation (49b) for L c is the downward velocity of the 
three-quarter-chord point of the airfoil. Since Up is the relative velocity 
conponent perpendicular to the quarter-chord, the sectional lift can also be 
written as 


%C 


1 c* c 

-pa - (-Up + - c) 

2 4 4 


(50a) 


L C 


1 c . 

- pacU (-Up + - e) 

2 2 


(50b) 


where V(t), appearing outside the parentheses of equation (49b), has been 
approximated by the resultant of only the tangential and perpendicular 
velocity components and is given by 


V a U “ U 2 + U 2 (51) 

T P 

As indicated in figure 5, the noncirculatory lift acts normal to the 
section chordline* and the circulatory lift acts normal to the resultant 
velocity U. The profile drag force acts parallel to U and is given by 


' *A portion of % c acts at the 3/4 -chord point and another at the 
1/2-chord point. However, the resultant of these two conponents is shown 
along the z 3 axis in figure 5 only for pictorial convenience. 
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D 


(52) 


1 Cd o ‘ 

D = - pac U 2 

2 a 

where is the (constant) profile drag coefficient* 

The conponents of the aerodynamic forces in the direction of the y 3 and 
Z 3 axes are given by 

(53a) 

(53b) 


F Y3 13 “Lc sin a - D cos a 


F Z 3 ~ cos ot + Ljjc "* D sin a 
where, from figure 5, 

sin a = Up/U (54) 

cos a = U<ji/U 

and U is given by equation (51). The aerodynamic force in the X 3 direction is 
given by F X 3 and is a profile drag force which is a function of the radial 
velocity component U R . Following usual practice, this force component is 
taken to be zero. 

Substituting equations (50), (52), and (54) into equation (53) leads to 


, 1 , c . C *o 

y 3 = - pac (U 1 2 U P £ U T U) 

2 P 2 a 


(55a) 


1 c.c. c 2 “ ^o 

f z 3 = - pac [-UpUm + - U T e UP + (-) e U P U] (55b) 

2 2 4 4 a 

The noncirculatory and circulatory moments in equation (49c) and (49d) can be 
written in terms of U<p, Up, U, and e and assume the form 


1 c 2 . . 3c 

Hun = “ “ pac (•) ( — Up “ Ue + ““ £ ) 

2 4 8 


(56a) 


1 c 2 . 

Mr. = - - pac (-) 2 Ue 

2 4 


(56b) 


from which the total pitching moment M X 3 is given by the sum of equations 
(56a) and (56b) as 
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( 57 ) 


1 c * . . J c 

«x 3 = - - pac (-) (Ue -Up + — e) 

J 2 4 8 

The virtual work of the aerodynamic forces can be written as 

6W A = j* t(F x e x + + F ? e 5 ) • 6R, 

0 

+ M X 3 [6<J) - (v' cos 0 + w* sin 0)(-5v' sin 0 + 6w' cos0)]]dx (58) 

where is the variation of the position vector of a point on the elastic 
axis and is given by 

SR-, = (6u - <$U F )e x 

+ (6v cos 0 + 6w sin + (-6v sin 0 + 6w cos 0)e^ (59) 

The components of the aerodynamic force vector F x , F^, and F^ are given by 



where [T] is the rotational transformation matrix which relates the coordinate 
axes of the deformed and undeformed blade and F X 3 has been set to zero* 
Substituting equation (59) into equation (58) , and integrating by parts, one 
obtains 


5w A = f ( A u fiu + A v 6v + A w 6w + A^fiijOdx + B A 
0 

where 

A u “ F x 

Ay = F n cos 6 - F ? sin 0 + (T A v' ) ' 

+ [- m X 3 (v' cos 0 sin 0 + w' sin^©)]' 


(61) 
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F_ sin 0 + Fr COS 0 + (TaW 1 )' 

+ [^X 3 (v' cos^8 + w 1 sin 0 cos 0)] 


(62) 



A<J> - M x 3 

t A “ f A u 
x 

j R R R 

B a = -T a v'6v I - T a W'6w + m X 3 (v' cos0 + w* sin0)sin0<5u 

1 0 0 0 

R 

- ^x 3 (v' cos 0 + w' sin 0 )cos 06w 

0 

In order to obtain the explicit expressions for the generalized aerodynamic 
forces, the quantities F y 3 , F z 3 , and M x 3 must be known in terms of the 
dependent variables u, v, w, and <j>, and the geometric pitch angle 0. This 
requires that U T , Up, and e first be obtained in terms of these variables. 
These egressions are developed in Appendix D, and are given by the equations 
(D16) and (D20). The expressions for U T and Up in equations (D16) include the 
induced velocity v^ which can be calculated using blade element and momentum 
theories. Equations (D16) and (D20) in combination with equations (51), (55), 
and (57) are sufficient to obtain the generalized aerodynamic forces from 
equations (60), (61), and (62). 

SUMMARY OF EQUATIONS 

In the previous sections, expressions for 6V, 6T, and 6W have been 
obtained. Substituting these expressions and their associated boundary terms 
into equation ( 1 ) , there results the egression of the form 

1 1 £ 

/ {/ [( )6u + ( )<Sv + ( )6w + ( )6<j>] dx + B} dt = 0 (63) 

to 0 
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For arbitrary admissible variations 6u, Sv, 5w, and 6$, the four expressions 
in parentheses roust vanish individually as oust the assembly boundary terms 
denoted by B. The first condition will yield the four governing nonlinear 
partial differential equations for u, v, w, and $ and the second condition 
gives the associated boundary conditions at the ends of the beam* Since the 
control angle 0 C is assumed to be given, the equation associated with the 
control will not appear. The governing equations of motion and boundary 
conditions are summarized belor. 

Extension: 

• • M •• 

• • 

mu - me( v 1 cos 0 + w* sin 0) — 2m ft (v - e<}> sin 0) 

- mft 2 ( x + u - Up - ev r cos 0 - ew 1 sin 0) + mft 2 Pp C (w + e sin 0 

+ e<f> cos0 ) - {EA[u' + k 2 4> 1 0 ’ - e A (v" + <f>w n ) cos 0 

A pt 

+ e A (<j>v" - w” )sin 0] + E^Au'}' = mg cos + A u (64a) 

Chor<±*ise bending* 

mv - me<f> sin 0 - 2m^0p C w - m ft 2 (v + e cos 0 - e<J) sin 0) 

- {me[ft 2 x(cos 0 - <j> sin 0) + 2ftv cos 0]} 

• • • • 

+ 2m£2(u - Up - ev 1 cos 0 - ew f sin 0) 

- (T c v')' “ [GJ<j>(v' 9 cos 0 sin 0 - w' ' cos 2 0 

+ v'Q* cos 20 + w'0' sin 20)]' + {EAe A u* (((> sin 0 - cos 0) 

pt pt 

- EB 2 <j> , 0 1 cos 0 + v* * [El^fsin 2 © + <j> sin 20) + EI^^(cos 2 0 

pt 

- <|) sin 20 ) ] + w 1 1 [ (EI^ — EI^Msin 0 cos 0 + <f> cos 20 ) ] 

+ GJ^ 1 (v* cos 0 sin 0 + w* sin 2 0) + cos 2 0 + 1^ sin^OJv 11 }* 1 

- (TgV 1 )* + [mge cos t|j(-cos 0 + <J> sin 0) + mge(cos 0 sin 

+ 0p C sin 0 sin ^Mv* cos 2 0 + w* sin 0 cos 0)]' = -mg sin \J> + A v (64b) 


35 



Flapwise bending: 


’ *• * 

• , 

nw + me<(> cos 0 + 2mfi0p C v - (Tqw' ) 

- {me [ft 2 x (sin© + <j> cos 0) + 2Qv sin 0]}' 

- sin 0 cos 0 + v' ' sin 2 0 + v'0 1 sin 20 

pt 

- 0' w' cos 20)]' + {-EAe^uM^ cos 0 + sin 0) 

pt 

- EB^'O'^ sin 0 + w' ' [EI^ cos 2 0 + EI^^ sin 2 0 

+ <)){EI ?i r - EI nn )sin20] + v” [(EI^ - EI nn )sin 0 cos 0 
+ <|>(El££ - EI^Jcos 20] - GJ4>'(w' cos 0 sin 0 + v' cos 2 0)}'' 

+ [E*(I 5 £ sin 2 0 + I nr| cos 2 0 )w' ' ] 1 ' - (T g w')' - (TcW')' 

+ tmge cos t("4> cos 0 + sin 0) + mge(cos 0 sin i]/ 

+ 0p C sin 0 cos sin 2 0 + v' sin 0 cos 0)]' = -mggp C cos i|» 

- m^PpoX + A w 
Torsion: 

• • 

mk^ <() + ) cos 2 0 + me [^^(w f cos 0 - v* sin 0) 

m ni2 m^j 

• • M 

- (v - ft 2 v) sin 0 + w cos 0] + mft 2 e<|> (vcos0 - x& pc sin0) 

- 2mft [e sin 0(u - U F ) - (k 2 - k 2 )v’ sin 0 cos 0 

m 2 m-| 

- w'(k 2 sin 2 0 + k 2 cos 2 0) - e0 DC (v cos 0 + w sin 0)] 

m 2 m-, * 

•• •• 

+ 2mftev(w' cos 0 - v' sin 0) - me<j>(v cos 0 + w sin 0) 

- {EAk 2 u' (0' + <j>' ) + EBi 0' 2 <f>' + EB 2 [0' ( <j>v* ' - w" ) sin 0 

& Pt pt pt 

- <j>'v" cos 0 - <j>'w' ' sin 0-0' (v* ' + 4>w' * )cos 0] + GJ r* • 

Pt 

+ (v'v' ' - w'w' ' ) cos 0 sin 0 - v'w' * cos 2 0 + w'v' ' sin 2 0 


(64c) 
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y 1 2 W* ^ 

+ 9' ( cos 20 - cos 20 + v’w' sin 20)]}' 

pt 2 2 

+EAe A u'(v' ' sin 0 - w" cos 0) + EB 2 (v"0' <j>' sin 0 - 

Pt 


- w' ' 0 ' <|> * cos 0 ) + v' 'w' ' (El rr ** EI^) cos 20 + w' ,2 (EI^^ 

pt 

. • 

- EI nn )sin 0 cos 0 + v" 2 (EI nn - EI cc )sin 0 cos 0 - [G J<|»’ 


+ E*Bi0' 2 4* * 1 * “ mge cos ♦(▼' sin 9 “ w' cos 0) 
pt 

- mge<J>(cos 0 sin + 0p C sin 0 cos ♦ ) “ mge(sin 0 sin i|> 

- 0p C cos 0 cos t|>) - mfl 2 0p C xe cos0 - mfl 2 sin0 cos 0 (k 2 - k 2 ) + (64d) 

The assembled collection of boundary terms denoted by B is given by 

B = Br - By + B$ W + B a + B a (65) 

D * 

and the requirement of the vanishing of the individual variational components 
leads to the relations 


( S 1 + d i “ m X 3 (v' cos 0 + w' sin 0) sin 0] 6u =0 

0 


(k ~ Trv' + s' — S 7 + s' + g - T a v' + d' - T A v' )6v “ 0 

2 C 4 ' 8 5 g 3 A lo 


( S 4 + S 8 + d 3) fi V 


[-k - Tpw’ + s' - s +s' + g - T a w' + d' - T A w 

4 5 9 10 6 g 4 


- M v (v' cos 0 + w' sin 0) cos 0] 6w 


£ 


( 66 ) 


( S 5 + s 10 + d 4)6w' 


( s 3 + s 6 + d 2 + d 5> 6 <t> - 0 

0 
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The tension T c and gravity force T g appearing in equations (64) and (66) are 
given by 

R ” •• 

T c = / m[”(u - Op) + e(v' cos 6 + w' sin 8) + 2ft(v cos 8 

x 

- sin 8) 

+ 8 2 (x + u - Up - ev' cos 8 - ew' sin 8) 

- 0 2 gp C (w + e<p cos8 + e sin8)] dx (67) 

Tg = mg cos ij> dx (68) 

x 

• • 

The terms Up and Up in the egressions for T c given in equation (67) lead to 
third-degree nonlinear terms when T c is substituted into equations (64) and 
(66) and can be discarded. Also, after substitution for T c in these 
equations, only resulting terms which are consistent with the ordering scheme 
adopted in Appendix C should be retained. Using the results given in 
equations (67) and (68) in combination with the extensional equation of motion 
given in equation (64a) (with dampings and A u set to zero), an alternative 
definition for the sum of T c and T g can be given as 

T C + T a = EA [u' + k 2 <|) ' 6 ' 
y A pt 

- e A (v' ' + <j>w' ' ) cos 8 + e A ( <}>v' ' - w" )sin 8] (69) 

The underlined terms in equation (64d) are associated with u' and they 
are called as tension-torsion coupling terras. These terras are known to be 
important in some cases (refs. 11 to 13). Also, there is a tension-bending 
term which is doubly underlined in the torsion equation. Some of these terms 
appear to be nonlinear, but they are not because of the relation given by 
equation (69). To simplify the solution of the equation (64), it is a 
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customary practice in the rotor blade literature to eliminate the extensional 
equation. If such practice is followed in the case of wind turbine blade 
solutions, one should substitute for EAu' in equations (64b), (64c), and (64d) 
in terms of T c + T g from equation (69). Also, after substitution only the 
resulting terras which are consistent with the ordering scheme adopted in 
Appendix C should be retained. 

As stated earlier, the generalized aerodynamic forces A u , A^, A w , and A<|> 
are obtained from equations (60), (61), and (62), using equations (D16) and 
(D20) in combination with equations (51), (55), and (56), and retaining the 
terras through second degree in the dependent variables. Because of the 
generality of the present development, these second-degree expressions are 
extremely lengthy and will not be shown. 

CONCLUDING REMARKS 

The second-degree nonlinear aeroelastic equations of motion for a 
flexible, twisted, nonuniform horizontal-axis rotor blade undergoing combined 
flapwise bending, chor&fise bending, torsion, and extension have been derived 
using Hamilton 1 s principle. The equations have their basis in the geometric 
nonlinear theory of elasticity and are consistent with the small deformation 
level of approximation in which elongations and shears (and hence strains) are 
negligible compared to unityi A mathematical ordering scheme which is 
consistent with the assumption of a slender beam was adopted for the purpose 
of systematically discarding elastic and <fynamic terms which are higher-order 
in the resultant equations of motion. The expressions for generalized 
aerodynamic forces, which account for windshear, tower shadow, and gusts, are 
left in general second-degree form from which one can obtain the aerodynamic 
forces loading to the order appropriate to any case of interest. A unique 
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feature of this development is the consideration of the pretwist of the blades 
before the elastic deformations which is more realistic than the common 
practice in the most of published work in which the pretwist is combined 
with the elastic twist. 
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APPENDIX A 

NONLINEAR CURVATURES AND COORDINATE SYSTEM TRANSFORMATION 
OF TWISTED ELASTIC WIND TURBINE BLADES, 

The nonlinear curvature ejqpressions for twisted rotor blades were 
developed in reference 2 in a general manner using the geometric nonlinear 
theory of elasticity. For convenience, this development followed a practice 
which is common in helicopter blade literature, namely the combination of the 
pre-twist with the elastic twist. As a consequence of enploying this 
simplification, if warping is considered, axial deformation exists in the 
initial configuration before any deformations are imposed. Such a situation 
would exist if an untwisted blade is twisted and then "frozen" to arrive at 
the pretwisted configuration. Since this situation does not arise in the 
fabrication of either wind turbine blades or helicopter blades, there is no 
axial deformation due to warping in the initial configuration before any 
deformations are inposed. In view of this situation, this Appendix will 
include the rotation due to pretwist with the control inputs and impose this 
rotation first while inposing the rotational transformations and will develop 
the second-degree nonlinear egressions for the curvature conponents and for 
the transformation matrix in terms of the variables u, v, w, and <|>. The 
foreshortening of the blades due to bending is explicitly considered. The 
level of approximation used within the geometric nonlinear theory of 
elasticity is that designated small deformations I in reference 2 in which the 
elongations and shear are negligible compared to unity and the square of the 
first derivative of the extensional deformation on the elastic axis is 
negligible coitpared to unity and the squares of the bending slopes. 

A schematic representation of the deformed and undeformed blade 
associated with a flap-lag-torsion transformation sequence is shown in figure 
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3 * The coordinate axis x is tangential to the elastic axis, and the axes t \ 
and £ are the principal axes of the section before deformation. Since the 
pretwist 0p t is not a constant, the initial curvature of the elastic axis 
before deformation is 
- d6pt 

w xns “ e x ~ e x®' 

ds pt 

where s is the distance measured along the undeformed elastic axis. 

The elastic deformations rotate the triad xnC to x 3 y 3 z 3/ and the 
transformation matrix between the two triads can be written as 

m 1 

m 2 n2 

m3 n3 

Let the expression for the curvature vector of the deformed elastic axis be 

^ X 3Y3 Z 3 = Wx 3 e *3 + W y3 e Y3 + Uz 3 ez 3 (A3) 

The next step is to find the expressions for the components of the 
curvature vector in terms of the direction cosines lj, m*|, • « • n3, and the 
blade initial curvature. From equation (A 2 ) one can write 

e X3 = l! e x + m-| e n + n-j e C (A 4 ) 

Differentiating equation (A 4 ) with respect to s and ejq>ressing the resulting 
expression along X3Y3Z3 axes, one can obtain 

e i = djl* + m^m* + n^n' ) e X3 

X3 11 1 

+ [1*l2 + m'm2 + n<n 2 ^ n 1 m 2 * m 1 n 2^ e Y3 

1 1 1 pt 
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+ [IT 3 + m'm 3 + n'lvj - 0 ' (n 11113 - 1111113 ) ] e Z 3 
11 1 pt 


(A5) 


Also the identity 


6v W x v Z ^ ®x 

X 3 3 y 3 3 x 3 


(A 6 ) 


leads to 


6y 3 Ww 6y “ Wy e z 

X 3 3 y 3 y 3 3 


(A7) 


From the identity 


Wy 67 ™ CO™ 6y 

X 3 Z 3 Z 3 X 3 


(A 8 ) 


one can write 


“x, " e * e z, 

3 y 3 3 


(A 8 ) 


Substituting equation (A2) into (A9) and using the orthogonal property of the 
matrix [T] , the expression for ui x ^ in terms of direction cosines is 


U) x = I 2 I 3 + n>2 m 3 + n 2 n 3 + ^1®' 

3 pt 

From equations (A5) and (A7) the expressions for w v and u> z are 

*3 3 


(A10) 


W, 


- [I1I3 + mim3 + nin3 - 0' (nim3 - 111^3)] 

pt 


u) z = [I1I2 + ®1 m 2 + n 1 n 2 “ ®' (»1“2 “ ®in2>] 


(All) 

(A12) 


"3 ' " ' _ ' " pt 

Thus for the expressions for the curvatures in terms of the direction 
cosines and the total section pitch of the blade have been developed* The 
next task is to express the direction cosines in terms of u, v, w, and <p . To 
this end, the direction cosines are first expressed in terms of the 
Eulerian-type angles $i, 5i, and 0i* It was shown in reference 2 that the 
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form of the expressions for the components of the curvature depends on the 
order in which the rotational transformations between the deformed and 
undeformed coordinates are imposed. In the present development, out of the 
six possible rotational transformation sequences which may be employed, a 
f lap-lag-torsion rotation sequence will be addressed. For this rotational 
transformation sequence, the rotations are imposed as follows: 

1. A positive rotation 3-j about the negative n axis resulting in 

x lYl z 1* 

2. A positive rotation C i about z^ axis resulting in X2y2Z2« 

3. A positive rotation Q-j about X 2 axis resulting in X 3 y 3 Z 3 . 

The eiqplicit form of the transformation matrix [T] in terms of the 
Eulerian-type angles 0i, Cl# and 0-) is 


CO - 


cos 0-| 
-sin Ci 
-sin 0i 
sin Cl 
-sin 0i 


cos Cl 
cos 0i 
sin 0i 
cos 0i 
cos 0i 


cos 0i 


sin 0i 


sin Cl 

cos Cl cos 0i 
-cos Cl sin 0i 


cos Cisin 0i 
cos 0i sin ©i 

-sin Cl sin 0i cos ©i 

sin Cl sin 0i sin ©i 

+cos 0i cos ©i 

(AI3) 


The rotation angles Cl# 01# 01 are to be expressed in terras of the 
variables u, v, w, <f>. To this end, let Rq and Ri be the position vectors of 
the point 0 and O' in figure 3, and AR be the displacement of 0. Then, one 
can write 


Ri = Rq + AR 

where 

AR = (u - Up) ex + vey + weg 

and 


(A14) 


(A15) 
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1 ,9 3v 2 3w 2 

U F - - / [(--) + (--) ] ds (A16) 

2 Q 3s 3s 


Using the following relation between the unit vector triads egeyez and ®x e ri e ? 


e X * e x 

ey = e^ cos 0 - e^ sin 0 (A17) 

e Z = e n 8 ^ n 9 + e x, 008 0 

the expression for AR is 


AR = (u - U F ) e x + (v cos 0 + w sin ©Je^ + (w cos 0 - v sin 0)e^ (A18) 

Differentiating equation (A18) with respect to s 

- 2 2 

3 ( AR) 3R oty a z 

"as"" = " + (a x “ " " " ) e x + a y e n + “z® 5 (A19) 

where 

°x " 

Oy = v' cos 0 + w' sin 0 (A20) 

a z * -v' sin 0 + w' cos 0 

1 s 2 2 

U F “ - / (oty + a z ) ds 

2 0 

Differentiating equation (A14) with respect to s and substituting equation 
(A19) into the resulting expression, the expression for 3R^/ds is 


- 2 2 
0 R 1 9 R Oty 0 Z _ 

'V. ■ Vs * <0 * ' 1 ■ 1 1 + “y*" + 

From calculus we have 


(A21) 
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(A22) 


3R 

3s 


ex 


Substituting equation (A22) into (A21) gives 


2 

3Rl CXy 

= (1 + 0 * 

3s 2 


2 

“z - 

"I ) e x + “yen + “z e S 


(A23) 


The relation between the extensional conponent of the Green's strain 
tensor e on the elastic axis and 3R-|/ds is given by 


3Rl 3R-| 

e = 1/2 ( . 1) (A24) 

3s 3s 


Substituting equation (A23) into (A24), the expression for e reduces to 



(A25) 


Invoking small strain assumption that the elongations and shears are small 
conpared to unity, equations (A23) and (A25) reduce to 

2 2 

3Rl Oy 0t z _ 

= (1 - — ) e x + Oye n + a z e? (A26) 

3s 2 2 


If the deformed length of the element ds is ds-j, the relation between d s«j 
and ds can be written as 

1/2 

ds 1 = (1 - 2e) ds (A28) 

Hence, 


3Rl -1/2 3R 1 

= (1 + 2e ) 

9s-| 3s 


Also, 


(A29) 
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(A30) 


3R 1 

ds-| 



Invoking small strain assumption and combining equations (A26), (A29), 

and (A30) results in 

2 2 

- “y a z - - 

e x ^ = (1 e x + a y e r) + a z e^ (A31) 


Prom equations (A2), (A13), and (A31) one can write 


ll = cos $i cos Si = 1 - 


2 2 
Oy 0 Z 

2 2 


mi = sin S i = Oy 

ni “ cos Si sin 0i = a z 
Hence, 


sin 0 i = a 2 


cos 0 i = 


1 - — 
2 


sin Si = dy cos Si = 


2 

, .i 

2 


(A3 2) 


(A33) 


The third rotation angle 0i is due to torsion of the blade and hence is given 
by 


®1 ■ ♦ . (A34) 

The egressions for the transformation matrix [T] and for the curvature 
components “X 3 , ^ 3 , u> z 3 are given in terms of the direction cosines in 
equations (A2), (A10), (All), and (A12) and those for the direction cosines in 
terms of the rotation angles are given in equation (A13). The rotation angles 
are ejq>ressed in terms of a x , cty, and a z in equations (A3 3) and (A34) and the 
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expressions for c^, oiy, and a z are given in terms of u, v, w, and <|> in 
equations (A20). Combining equations (A10), (All), (A13), (A33), (A34), and 
(A20), the second-degree expressions for the curvature components and for the 
transformation matrix can be obtained as 


u x 3 


<j>' - (v' cos 0 + w' sin 6)(-v' sin 0 + w' cos 0)' 


v' 2 w' 2 

+ 0 ' (1 ) 

pt 2 2 


(A3 5) 


u y 3 = -w' ' (cos 0 - <j> sin 0) + v' ' (sin 0 + <J> cos 0) (A36) 

a) Z 3 = v* ' (cos 0 - <J> sin 0) + w' ' (sin 0 + $ cos 0) (A37) 



v' 2 w' 2 

1 v 1 cos 0 + w* sin 0 -v' sin 0 + w* cos 0 

2 2 


-v* (cos 0 - <J> sin 0) 
-w*(sin 0 + <J> cos 0) 


1 - 


(v 1 cos 0 


+ w 1 
2 


sin 0) 


<f > 2 

2 


<j> - (V* cos 0 + W 1 ) 

• (-v 1 sin 0 + w 1 cos 0) 


v* (sin 0 + <j> cos 0) -<|> 

-w' (cos 0 - <|> sin 0) 


1 


(V 


(j ) 2 


sin 0 + w* cos 


2 


2 


(A3 8) 
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APPENDIX B 


NONLINEAR STRAIN-DISPLACEMENT RELATIONS 
This Appendix will develop second-degree nonlinear expressions for 
strains* To this end, let Tq and r^ be the position vectors before and after 
deformation of an arbitrary mass point on the blade* These vectors can be 
expressed as 


r 0 = Rq + ne n + Ce^ 


(Bl ) 


3^1 ~ R^ts^) + ne y + 

**3 3 


(B2) 


where s-| is the length measured along the deformed elastic axis of the blade* 
The differentials of the above vectors are given by 


drQ = dse x + (dn “ S© 1 ds) e^ + (d^ + n9' ds) e^ 

pt pt 


(B3 ) 


dr-j = (1 - e x ^) ds-j 

+ (dn - ds-|) e y ^ + (d£ + no) x ^ ds-|) e z ^ (B4) 

where the curvature corrponents u> x , o> y , w z are defined in Appendix B* The 

3 3 3 

usual practice in solid mechanics is to use the Lagrangian strain tensor which 
is defined by 


ds 

dr 1 . dr-j - dr 0 . dr 0 ■ 2 [ds dn d£] [e^j] ^ dn 

d£. 


(B5) 


Substituting equations (B3), (B4), (A28) into equation (B5) and using the 
relation between the engineering strains and the components of the Lagrangian 
strain given by equation (5), one obtains the following expressions for the 
strain components 
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<J) ’ 2 

Yxx " e xx = u ' + (n 2 + 5 2 ) (--- + ♦' 8^) 

- n[v' ' cos 0 + w' ' sin 0 + <j) (-v‘ ' sin 0 + w' ' cos 0)] 

+ 5 [v' ' sin 0 - w' ' cos 0 + <J> (v 1 ' cos 0 + w' ' sin 0 ) ] (B 6 ) 


Y xn = 2e xr) = -C t<J>' " (v' cos 0 + w' sin 0) 


v ' 2 w '2 

• (-v' sin 0 + w' cos 0)' -0' ( + — )] (B7) 

pt 2 2 


Y X £ = 2 c x ^ = n t 4 >’ - (v' cos 0 + w' sin 0 ) 


(-V 1 sin 0 + w' cos 0)' -0' 

pt 


v »2 w >2 



(B8) 


It should be pointed out that in arriving at the expression given in 
equations (B 6 ), (B7), (B 8 ) several higher-order terms have been discarded 
based either on considerations related to small deformations in which 
elongations and shears are small conpared to unity or on considerations 
related to the approximations which can be made because of the slenderness of 
the blade as discussed in Appendix C. Retention of higher order terms in the 
expressions for the strain conponents is not a problem. However, these higher 
order terms in the strains lead to higher order terms in the final equations 
of motion. Thus, discarding these terms in the resultant strain expressions 
using the considerations of Appedix C sinplifies the algebric manipulations. 
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APPENDIX C 


SLENDER BEAM APPROXIMATION AND ATTENDANT ORDERING SCHEME 
The simplifications of the slender beam approximation as applied to the 
derivation of the second-degree nonlinear equations of motion were discussed 
and a mathematical ordering scheme which is compatible with the assumption of 
a slender beam was introduced in reference 3* The same scheme has been used 
in the present report* In this scheme, a slender beam is systematized by 
introducing a parameter e which is taken to be of the same order as the 
nondimensionalized bending displacements v/R and w/R. The order of the 
dependent variables and the geometric quantities appearing in the equations of 
motion of this report are as follows: 


u/R 

- 0(e 2 ) 

n/R = 

0(e) 

e c " 

0(1) 

v/R 

= 0(e) 

S/R = 

0(e) 

6' = 

0(e) 





pt 


w/R 

= 0(e) 

O 

1! 

0(e) 



<f> = 

0(e) 

x/R = 

0(1) 





9 pt - 

0(1) 




Following this ordering scheme, the order of the elastic and inertial 
terms which are retained in the second-degree nonlinear coupled 
flap-lag-axial-torsion equations of motion of this report are given in Table 
Al below* 

TABLE Al. - ORDERING SCHEME 
Freedom Elastic forces Inertial forces 

Bending e 4 e 2 

Torsion e3 

Extension e3 e 3 

The rationale for this scheme was discussed in reference 3. 
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APPENDIX D 


VELOCITY COMPONENTS EXPERIENCED BY BLADE ELEMENT 
The velocity conponents experienced by a blade element are determined by 
considering the contributions from the atmospheric winds, induced velocity, 
and blade (fynamic velocity. For a down wind rotor, the tower shadow effect 
will be included. 

The resultant velocities seen by a point on the elastic axis of the blade 
in deformed and the undeformed coordinate systems are related according to 

’ V! .J - m ''xnc <D1 > 

where, from figure 6, 


Vx 3 y 3 Z 3 = Ur6x 3 ~ UTey 3 ~ UpSz 3 


(D2) 


dr! 

and V xn? = (V a - ~-) xn5 (D3) 

at 

The aerodynamic velocity conponents seen by a blade element are shewn in 
figure 7. The aerocfynamic velocity consists of three conponents: (1) the 

free stream velocity; (2) the gust velocity; and (3) the induced velocity. 

The free stream velocity profile over a rough terrain is frequently 
approximated by a power-law relation with height, and is given by 

xj n 

V m (x^ = V H (1 - --) (D4) 

H 

where V m and V H are the mean velocities at xj and turbine axis respectively 
and H is the height of the center of the hub from ground. This mean velocity, 
V m is a function of the azimuth angle of the blade since the interference of 
the tower causes a reduction in the velocity when the blade is in the 
vicinity of the tower. In the present development, it is accounted for by 
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introducing a blockage factor b which is less than or equal to 1. Inside the 
tower shadow band defined by the limits ^ < ^2 the roean velocity is bV m 

where is defined in equation (D4). The limits i)>i and <Jj 2 and the blockage 
factor b are to be obtained from either wind tunnel experiments or full scale 
tests* Expanding V(xj) in Taylor's series about xj * 0 and retaining only up 
to second-order terms/ equation (D4) reduces to 

“ X H (1 + Ci x 2 “ C 2 X cos i|» + Ci X 2 cos 2i|>) (D5) 

where 

V - V m /flR 

Xj = X COS $p C COS “ X COS l|» 

Ci = (n 2 - n)R 2 /4H 2 

C 2 = nR/H (D6) 

X =■ x/R 
X H - V H /flR 


The gust velocities are assumed to be harmonic with a frequency U)g and 
are represented as 


u g 

Xug * — sin o) g t 

JlR 


v g 

X^g “ sin u)gt (D7 ) 

OR 


Xwg = — Sin OJgt 
SiR 


The velocity V m = X m ftR is inclined to the rotor axis at an angle as 
shown in figure 7. Resolving this velocity along Yj and Zj axis, one can 
write 
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(D8) 


v m = bv m sin + bV m cos ^oez^. 

= ftR (bX m sin ipoey^ + bX m cos \|>oe z ^) 

The induced velocity can be calculated using blade element and momentum 
theories and is 

Vi = XiflR (D9) 

How the aerodynamic velocity can be expressed as 


VX I Y I Z I = ^ CXu 9 ex I + (bXm Sin + X vg)eY I 

+ (bX m cos \J»o + X wg - XiJez^] (DIO) 

This velocity can be expressed in xr)£ system by the relation 

v xnC = (Dll) 


where the transformation matrix can be written with the aid of the figures 1 
and 2 as 


PJ 


COS 

“®pc cos ^ sin 0 
-sin cos 0 

-3 PC cos cos 0 


+sin sin 0 


sin ij# 

“frpc sin sin 0 

+ cos ij; cos 0 

*$pc sin cos 0 

-cos sin 0 


Ppc 

sin 0 


cos 0 


(D 12) 


The position vector of a point on the elastic axis is 
r-j = (x + u - Up)e x 

+ (v cos 0 + w sin 0)e n + (-v sin 0 + w cos Q)e^ (D13) 

The angular velocity of the coordinate system xr}£ is given by equation (20), 
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Using equation (D13) in conjunction with equation (20) , the dynamic velocity 
of a point on the elastic axis taken with respect to the xn? axis system is 
given by 
dri 

= (u - Up - ftv)e x 

dt 

+ [v cos 0 + w sin 9) - ft£p C (-v sin 0 + w cos 0) 

+ ft cos 0 (x + u - UpJJe^ 

+ [-v sin 0 + w cos 0 + &3p C (v cos 0 + w sin 0) 

- ft sin 0 (x + u - Up)]e£ (D14) 

Combining equations (Dll), (D12), and (D14) gives the total velocity seen by a 
point on the elastic axis as 

v xns = cos t + (bX m sin i|>0 + X vg )sin t 

+ 0 pc (bX m cos to + Kg ~ *i)l 

- u + Up + 8v}e x 

+ cos t cos 0 - sin t sin 0) 

+ (bX m sin to + X vg )(-0p C sin t sin 0 + cos t cos 0) 

+ (bX m cos 4»o + *wg " X^Jsin 0] 

- (v cos 0 + w sin 0) + Q0p C (-v sin 0 + w cos 0) 

- Q cos 0 (x + u - Up)}e T1 

+ {fiR[X U g(-0p C cos t cos 0 + sin t sin 0) 

+ (bXj,, sin to + X V g)(-0p C sin t cos 0 - cos t sin 0) 

+ (bX m cos to + X wg - Xi)cos 0] 

• • 

- (-v sin 0 + w cos 0) - ^3p C (v cos 0 + w sin 0) 

+ ft sin 0 (x + u - U F )}e^ (D15) 

The tangential and perpendicular velocity components Uij» and Up are 
obtained from equations (A38), (Dl), (D2), and (D15) and, to second-degree in 
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the dependent variables, have the form 

U T = - {-RRlA^jg cos + (bX m sin iJjq + X V g) sin ij) 

Ppc cos +0 ^ \/g ” X^)] 

. [v' (cos 0 - <f> sin 0) + w' (sin 0 + <J> cos 0)] 

+ (-u + ftv) (-v* cos 0 - w' sin 0) 

+ fiR[X ug (-3p C cos ij; cos 0 - sin sin 0) 

+ (bX m sin i|»o + X V g)(-0p C sin ip sin 0 + cos )Jj cos 0) 

+ (bXfl, cos iJiq + X W g - X^) sin 0] 

(v' cos 0 + w' sin 0)2 <b 2 

. M :-] 

2 2 


- (v cos 0 + w sin 0) + fiPpc (-v sin 0 + w cos 0) 


. . , (v' cos 0 + w' sin 0) 2 <fc 2 

- ft cos 0 (u - U F ) - ftx cos 0 [1 — ] 

2 2 


+ ftR [X U g(-3 pc cos ij> cos 0 + sin ij> sin 0) 

+ (bX m sin \J>o + X vg )(-3 pc sin ♦ cos 0 - cos ^ sin 0) 

+ (bX m cos + X W g - Xj_ ) cos 0] 

• [<f> - (v' cos 0 + w' sin 0)(-v' sin 0 + w' cos 0)] 

• • 

+ ♦ [” (“V sin 0 + w cos 0 - ft3 pc (v cos 0 + w sin 0) 

+ ft sin 0 (x + u)]} (D16a) 


Up = - {ftR [X U g cos i|> + (bX m sin i[>q + ^vg) sin ♦ 

+ 3 pc (bX m cos i|) 0 + Xv, g - Xi>] 

. [v' (sin 0 + <J> cos 0) - w' (cos 0 - <|> sin 0)] 

+ (-u + ftv) (v f sin 0 - w' cos 0)] 

- ftR [Xug (-3 pc cos ip cos 0 - sin sin 0) 

+ (bX m sin i|»o + X V g) (~3p C sin sin 0 + cos cos 0) 

+ (bXjj cos \p 0 + *wg _ *i> sin 0] <(> 

- [-(v cos 0 + w*sin 0) + ft3p C (-v sin 0 + w cos 0)- ft cos0 (x+u)]<J> 
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+ ftR [\ig(“3pc cos ^ cos 0 + sin i|» sin 0) 

+ (bV sin ij»o + A vg )(-e pc sin ♦ 008 9 " 008 + sin 9 > 

+ (bX m cos \|»o + " ^i) cos 0] 

- (-v sin 0 + w cos 0) - S0p C (v cos 0 + w sin 0) + ft sin 0 (u - Dp) 

(v' sin 0 + w' cos 0)2 <J> 

+ !) sin 0 x (1 )} (Dl6b) 

2 2 

• 

The quantity e appearing in equations ( 51 ) and ( 53 ) is the angular 
velocity of the blade section about the local X3 axis and, consistent with the 
present notation, can be written as e x ^. It can be regarded as composed of 
two parts: the first part arising from the rigid-body angular velocity of the 

hub in space, the second part arising from the angular velocity associated 
with the elastic deformations. The first contribution due to 0 and is given 
by 


Then 


^ 


r 

% 

i j 

O 

c; 

< %\ 

^ = (T] > 

1 ! 
ft sin 0 

< e z - 


^.ft cos 0^ 

*5 

ft 



v' 2 

w *2 

S = 

fiPpc (1 - --- - 

) + ftw 

2 


(D 17 ) 


(D 18 ) 


The second part of the angular velocity of blade section can be obtained by 

• I • 

replacing <(>' by <j> , (-v' sin 0 + w' cos 0) by (-v' sin 0 + w' cos 0), and 
0'pt by 0pt (which is zero) in the expression for the curvature io x ^ given in 
equation (A 35 ). Thus, 


(e x ^) deformation = 5> " (v* cos 0 + w' sin 0 )(-v' sin 0 + w' cos 0 ) (D 19 ) 
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n • 


Combining equations (D18) and (D19), the total sectional pitching velocity is 


v'2 w'2 

x 3 > - a 6pc - --> ♦ aw' 

+ - (v 1 cos 0 + w 1 sin GM-v* sin 0 + w 1 cos 0) 


(D20) 
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Figure 2. - Coordinate systems of blade cross 
section. 


gure 3. - Schem; 
deformed and de 
igle, 0, not shoi 



/ 



65 


Figure 4. - Cross section of blade in 
general unsteady motion. 
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Figure 5. - Blade section inflow geometry and 
aerodynamic force components. 
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